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Abstract 

In this work, we study the charged lepton flavor violating (cLFV) decays r — > fij, 
t — > 3/i and Z — >■ jit in the framework of the supersymmetric economical 3-3-1 
model. Analytic formulas for branching ratios (BR) of these decays are presented. 
We assume that there exits lepton flavor violation (LFV) sources in both right- and 
left- handed slepton sectors. This leads to the strong enhancement of cLFV decay 
rates. We also show that the effects of the LFV source to the cLFV decay rates in 
the left-handed slepton sector are greater than those in the right- handed slepton 
sector. By numerical investigation, we show that the model under consideration 
contains the relative light mass spectrum of sleptons which satisfies the current 
experimental bounds on LFV processes in the limit of small tan 7. The interplay 
between monopole and dipole operators also was studied. 

Key words: Decays of taus, Decays of Z bosons, Supersymmetric models 
PACS: 13.35.Dx, 13.38.Dg, 12.60.Jv 



1 Introduction. 



It is known that lepton flavor (LF) number are strictly conversed in Stan- 
dard Model (SM). However, the observation from of neutrino experiments [5] 
strongly suggests that there is lepton flavor mixing in lepton sector. It means 
that charged lepton flavor violating processes may also occur at some level. 
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Many current experiments especially pay attention to search for LFV pro- 
cesses in the charged lepton sector such as tau decays [1,2], Z — \ \it decays 
[3,4], ect... because LFV observed in experiment is evidence of new physics be- 
yond the SM. Along with experiments, many models beyond SM that contain 
LFV processes have been constructed. One class of the simple extended SM 
models with LFV is the class of models with non-zero neutrino masses. This 
kind of models contains a new type of Yukawa couplings of right-handed neu- 
trinos which are sources of cLFV. But cLFV processes in these models have 
been proved to be very suppressed [21]. However, in supersymmetric models 
(SUSY), the situation can be changed. Besides LFV origin affected by the new 
Yukawa couplings involving right-handed neutrinos, SUSY models also con- 
tain the other source of cLFV. Particularly, the large mixings of slepton mass 
parameters in the soft term greatly enhance the rates for cLFV processes. 
Because of the appearance of new cLFV sources in SUSY models, the cLFV 
decay rates can be reach to experimental bounds even if the mixing angles in 
the slepton sector are small. Other interesting properties supporting the study 
cLFV phenomenology in SUSY models were discussed in many recent works, 
for example [26]. Many investigations about the cLFV in SUSY models such 
as [9,22,8] indicated that the values of BR(r — > 3/x) and BR(r — > /ry) may 
exceed the bounds of present experimental results. So it is necessary to find 
regions of parameter space satisfying the experimental results [2,3,4]: 

BR(t~ ->■ // _ 7) < 4.4 x 1(T 8 , (1) 

BR(r~ ->■ /x~a* + aO < 2.1 x 1(T 8 , (2) 

BR(Z ->■ fi + r-) < 1.2 x HT 5 . (3) 



Supersymmetric unified theories predict large rates for cLFV processes except 
the decay rate of Z — > fir. The branching of this decay in models beyond SM 
is predicted to 10~ 9 — 10~ 8 level [7]. It is very suppressed compared to the 
experimental bound (3). In fact, the GigaZ option will be approaching to the 
sensitivity of 10~ 8 level in the Z — > fir decay mode [6]. If observed in a future 
experiment, it will be evidence of physics in SUSY models. One of the SUSY 
models in which cLFV processes are thoroughly investigated is the minimal 
supersymmetric standard model (MSSM) [9,10]. In the mentioned work, the 
authors have showed that if there exists LFV in the left-handed sector, in 
order to get the experimental bound on the LFV decay rates of muon and 
tauon decays, the mass parameters of sleptons should be shifted toward TeV 
scale, especially in the case of large tan /3. Because the model contains many 
LFV effective operators, the interplay of different effective operators (dipole, 
monopole coming from neutral boson exchanges and Higgs exchanges) creates 
many interesting consequences. The detail is discussed in works [9,10]. 
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The current experimental results [5] show that neutrinos are massive, which 
contradict what assumed in the SM. Other words speaking, the SM must be 
extended. Among extensions, the models based on the SU(3)c®SU(3)l®U(1)x 
(3-3-1) gauge group [12,13] have the following interesting features: 

• To be anomaly free, the number of triplets should be equal to number of 
antitriplets. This leads to that number of generations is multiple of the color 
number (3). Combining with condition of the QCD asymptotic freedom the 
requiring number of quark generations should be less than five. Thus, in the 
the 3-3-1 models, number of generations is three. 

• The models give an explanation of electric charge quantization [14], dark 
matter and CP violation. 

• One of the three quark families has to transform under SU(3)l differently 
from the other two. This leads to an explanation why the top quark is 
uncharacteristically heavy. 

The weakness of the 3-3-1 models is the complication in the Higgs sector: in 
the minimal 3-3-1 model [12], content of Higgs sector includes three triplets 
and one sextet, while in the 3-3-1 model with right-handed neutrino [13], there 
are three triplets. To solve the mentioned weakness, the attempt is reduction 
in the scalar sector. The first result is the economical 3-3-1 model [15,16] - the 
version with right-handed neutrino and only two Higgs triplets. Scalar sector 
in the minimal 3-3-1 model has been reduced to minimum with two Higgs 
triplets, and such version is called the reduced minimal 3-3-1 model [17]. It is 
emphasized that the problem on neutrino masses is not totally solved in the 
last version. The situation seems better in a supersymmetric version of the re- 
duced minimal 3-3-1 model [18]. An supersymmetric version of the economical 
3-3-1 model has been established in [19] and called the Supersymmetric Eco- 
nomical 3-3-1 model (SUSYE331). The SUSYE331 has the same advantages as 
its non-supersymmetric version. However, to generate masses for fermions, the 
non-supersymmetric models with minimal Higgs content need effective non- 
renormalizable interactions, while the SUSY versions do not (the interested 
readers can find in Refs. [19,24]). 

In [20], the LFV decays of neutral Higgs bosons in the SUSYE331 were consid- 
ered. In this work, we are interested in the cLFV processes in the SUSYE331, 
namely: cLFV decays of the tauon and Z bosons. 

We remind that in the SUSYE331, there are more leptons, Higgses and gauge 
bosons as well as their supersymmetric partners than those of the MSSM. 
This implies that the model contains many sources of LFV in the slepton sec- 
tor. This suggests that the branching ratios of LFV Higgs and Z decays are 
greatly enhanced and may be detectable in near future. In contrast, bounds 
from experiments (such as shown in (1) and (2)) will strictly constrains the 
values of parameters causing the cLFV. So it is really important to investi- 
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gate the parameter space where LVF decay rates can be satisfied the present 
bound of experiments. In many previous works [9,10], the authors showed 
that in the MSSM the interested region of parameter space is set at TeV scale 
with the large value of tan (5. Because of appearance of new LFV sources in 
the SUSYE331, we can find the interested region of parameter space even in 
the limit of small slepton mass parameters. In this work, we are interested 
in the processes: t~ — > /i~7, r~ — > fi~fi + fi~ and Z — > fir. In particular, we 
incorporate the mixing of sleptons mass, the mixing of charginos, Higgsino, 
gauginos as well as the interactions between gauge boson, slepton, lepton and 
the Yukawa interactions. This leads to many types of enhanced diagrams. We 
will consider how large contribution from each type of diagrams can be. We 
also will extend previous work in [9] to our considered model such as: Con- 
structing the analytical formulas of effective operators from the diagrams, from 
which we obtain the formulas of the branching ratios of decays BR(r — > /ry), 
BR(r — > fififi) and BR(Z — > fir). After that, we investigate some numerical 
results of these branching ratios in the limit of small tan 7 and the slepton 
parameters are set below 1 TeV which may be detected by current colliders. 
Here, the quantity tan 7 in the SUSYE331 plays a similar role as the tan/3 
in the MSSM [20]. In the MSSM model, the mass of the lightest Higgs de- 
pends on tan j3 and the mass of the standard gauge bosons. If we combine the 
theoretical result for the upper bound on the lightest Higgs-boson mass with 
the direct experimental search, it leads to exclude the limit of small tan (3 for 
the case where soft parameters are set to TeV scale. However, this kind of 
constraints on the tan 7 does not happen in the SUSYE331 model. On the 
other hand, the large values of tan 7 do not support the region of the soft pa- 
rameter space below 1 TeV. Hence, in this work, we concentrate on numerical 
studying in the limit of small tan /3. In this paper we concentrate on only two 
aims. First, we establish analytic formulas to calculating some cLFV processes 
in the framework of SUSYE331. Second, we prove that there exist regions sat- 
isfying the current bounds of cLFV experiments. Especially, on the basis of 
numerical studying, we find some interested regions of parameter space that 
satisfy the experimental bound on the cLFV decay rates. We also discuss on 
the interplay between mono and dipole operators which may lead to some 
interesting consequences. 



This paper is arranged as follows: In section 2 we review the particles content. 
The effective operators as well as branching ratios of the mentioned cLFV 
decay processes are presented in section 3. Section 4 is devoted for discussion 
on results of numerical investigation on parameter space and the last section 
is the conclusion. All the analytic formulas of effective couplings appearing in 
effective operators and interacting terms needed for our calculation are shown 
in Appendices. 
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2 Particles content 



In this part, let us quickly review the particle content in the SUSYE331 model, 
which were given in the previous papers [19,24,20]. The fermion superfields are 
given by 

L aL = (v a ,T a ,Vi) T L ~ (1,3,-1/3), 1^ - (1,1,1), a = 1,2, 3 (4) 

Qil= («i, d u u') T L -(3,3,1/3), (5) 

u c 1L , ££-(3^,-2/3),^- (3*, 1,1/3), (6) 

Q a L = (d a ,-u a ,d' a ) T L ~(3,3*,0), a = 2,3, (7) 

u c aL ~ (3*, 1, -2/3) , d c aL , fa ~ (3*, 1, 1/3) , (8) 

and the Higgs superfields are written as 

X = (Xi, r , X°) T ~ (1, 3, -1/3), p = (pf, fP, p+f ~ (1, 3, 2/3), (9) 
X = (x5r^ + ,X?) T ~ (1,3*, 1/3), ft = (ftr,ft°,fr) T ~ (l,3*,-2/3).(10) 

It is noted that ip c L = {ipR) c = ^j? an d are exotic quarks. The values in 

each parenthesis show the quantum numbers of the (<S77(3) C , SU(3)l, U(1)x) 
group, respectively. The SU(3)i ® U(l)x gauge group is broken as follows: 

SU(3) L ® U(l) x ^ SU(2) L ® U(l) y ^ U(1) Q , (11) 

with VEVs given by 

V2( X f = (u, 0, w) , V2( X 'f = (u\ 0, w') , (12) 
v / 2(p) T = (0, v, 0) , V2(p') T = (0, 0) . (13) 

The SUSYE331 model contains seventeen vector superfields such as V®, V a 
and V'. The vector superfields contain the usual gauge bosons given in [19,24]. 
The total supersymmetric Lagrangian was given in [24]. In this work only 
terms relevant to our work are collected in appendix A. 
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3 Effective operators and branching ratios 

In this section, we extend the previous work [9] in the model under considera- 
tion. In the SUSYE331, there are six physical vector bosons: the photon, two 
charged (W ± , Y ± ), two hermitian neutral (Z,Z') and a non-hermitian neu- 
tral X. The effective couplings of Xfir interaction is very small so we ignore 
them in the calculation. Let us first write down the cLFV effective operators 
for muon-tau and photon or Z bosons, Z' bosons or other leptons. 



3.1 77/7 effective operators 

First we write down the LFV operators for r, /x, 7. These operators are divided 
into two terms 0: 

e [Cljia^r + C>V M f c + h.c] nA„, (14) 
e m T [DlJla^f + D^a^r + h.c] F^. (15) 

The processes associated with external photon line depend only on the D\ R 
while other processes with virtual photon depend on both C 7 and D 1 . The 
Feynman diagrams contributing to C 7 and D 1 are given in appendix B. We 
would like to emphasize that the number of diagrams in the considered model 
is more than that of the MSSM because the SUSYE331 model contains new 
Higgs and gauge bosons. In order to obtain analytical formulas for C 7 and 
-D 7 , we have used Feynman rules and some approximate expansion. For more 
details, the interested reader can see in the appendix B. As in the MSSM, 
the C2(i?) does n °t depend on tan 7 and the factor D^rm can be divided into 
three sub-terms including sub-term DffijL which is independent on tan 7 and 
sub-terms -D^fe) which are not, more specifically 

Dl iR) = Dlfl+Dl^ )+ Dl^ y (16) 



We would like to mention that only D£^L contains left-right slepton mixing 
parameters which, in this paper, are denoted by (A T , A^ T , A^ T ). 



4 The operator containing D 7 is different from [9] a factor i. This is because of 
definition of and a^ v ' . The definitions in [11] are used. 
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3.2 Zrp and Z'rp effective operator 



First, let us consider the effective operator of the Zrp. This kind of effective 
operators can be written in the standard form given in [9] as follows: 



9z m z 



Alfia^T + Altfo»f c + h. 



(17) 



9z 



9z m T 



Clfiu^T + Cf/iV^f c + h.c 



UZ 



it- 



Dlfta^r + Dltfo^T + h. 



rP 



j fii/- 



(18) 
(19) 



where the mass of the Z boson in the SUSYE331 is determined as [19]: 



2 g 2 (v 2 + v' 2 ) 



Ac 2 



g z ^— and gz' — — - — • (20) 



{gz and g z > are defined from covariant derivatives, as explained in the Ap- 
pendix). The operators related to the factors Af R are chirality conserving 
-monopole. The leading contribution to the monopole operators comes from 
effective couplings of muon and tauon with two neutral Higgs bosons. In the 
model under consideration, we have four neutral Higgs bosons p°, p'°, Xi an d 
xf which can couple to the Z boson. However, investigation in [16] has noted 
the relation u, u' w, w' . In this limit, we can neglect the coupling of 

the Z with \ an d x' ■ It means that we obtain only the leading interactions of 
fj,,T with p,p'. These leading terms lead to a consequence that the monopole 
operators in (17) can be extracted out a factor m 2 z . We would like to remind 
that the class of diagrams containing /i, r, p° and p'° presented in appendix 
(C.l) give contribution to Af R and this factor can be written as sum of three 
parts: 

aZ _ A Z(a) ,Z(b) ,Z(c) / 91 x 

L(R) — ^UR) + A L(R) + L(R) ' \ ZL ) 



where analytical formulas of each term in (21) as well as Ct{R) an d ^l(r) are 
given in appendix C. The operator related to CftR) a I so are chirality conserving 
-monopole while the operators related to Cf(m a re chirality flipping (dipole). 
More details about origin of these operators, the interested reader can see in 
[9]. 
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A big difference compared to the MSSM is that the model under consideration 
contains more pr effective operators such as prZ' . Because of the couplings 
of Z' with all of neutral Higgs -X%x' 2 ^P an d p'- the monopole operators of 
Z'pr relate not only to the factor m| in but also the factor m|,, as indicated 
in two formulas (22) and (23). It is noted that A lZ ^ comes from the leading 
interactions of pr with two neutral Higgs p and p' while the A 1Z ^ comes 
from the leading interactions of pr with two neutral Higgs \\ an d X2 ■ For a ^ 
of others effective operators mentioned in this work-Cf'^ and D z '[ R y relate 
with only two neutral Higgs p and p'. Above comments are enough for us to 
write the standard form of prZ' operators as follows: 



g z , m 2 z [A^'Wr + Ag z '>/2 W + h.c] Z^ 



(22) 



9z> m z 



9z> 



, [Af'Wr 



+ A 



{ l z,) p c cj»f c + h.c 



Z' 



Cfpa^T + C z R p c a»f c + h.cj UZ'p 



(23) 
(24) 



9z' m T 



Df'pa^T c + D z R p c a» u T + h.c 



(25) 



Note that in the on-shell condition we have O(Z) — > —m 2 z and O(Z') — y - in 
where: 



2 

Z'i 



m|, 



An 2 



W 2 = w 2 + w' 2 . 



The AW r -> and are also written in the form as those for Af, R y Forms of 

the -Df(R) and Df'^ are the same of -DLm ^ n (^)- ^ analytical formulas of 
effective operators in this section are given in appendix D. 



In addition, the SUSYE331 model has two others neutral gauge bosons, Z 1 and 
X compared to the MSSM. This appearance may give significant contribution 
to the r — >■ 3p decay. The mentioned contribution may be similar to that of 
the Z boson in some region of parameter space. This issue will be discussed 
in more details in the next section. 



3.3 rppp effective operator 



Let us write down four-fermions rppp effective operators. The rppp effective 
operators can be constructed from the effective operators of p, r with photon 
or Z, Z' bosons as well as Higgs bosons. Besides the contributions coming from 
photon, Z boson and Higgs exchanges, there are other contributions to the 
rppp effective operators which come from box-diagrams shown in fig. (E.l). 
In this part, we write down only the standard form of rppp effective operator 
coming from box-diagrams as 
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+ (/i W) (B% + B£ fi c a,fi c ) + h.c] . (26) 



Analytical form of the coefficients B^j^ are presented in the appendix E. 



3.4 Branching ratios 



General method to construct the branching ratios for the cLFV decay from 
effective operator was written in [9]. Especially, based on the basis of the 
effective operators, we write effective Lagrangian of the muon and tauon with 
photon, gauge bosons Z, Z' as well as lepton. From this effective Lagrangian, 
we obtain the branching ratio for each of processes by using the Feynman 
rules. In this section, we study the branching ratios for the considered cLFV 
decays of the tau and the Z bosons. 



3.4.I t — y fj,^f 

Let us consider the cLFV decay mode r — > /ry. It is not hard to obtain the 
decay rate of r — >■ /ry from effective operators given in Eqs. (14,15). The 
detailed calculation can be seen in [22]. Comparing the branching of LFV 
decay r — > /ry with that of the r — > fiu^ decay given in [4,11], we can obtain 
the relation between the two branching ratios of the two above processes. The 
result is entirely consistent with the results given in [9], namely 



BR(r ->■ ^-7) 



48vr 2 a 



Dl\ 2 + \Dl\ 2 BR(r^^V T ), 



yy 1 2 



(27) 



where a = g, = ^- and Br(r -> ^v^) ~ 17.41%. 
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3.4.2 Z ->• /ir 

In this subsection, we consider the decay mode of Z — > fir. The decay rate 
of Z — > £ + £~ in case of the SUSY can also be determined from the general 
formula established in [11]. In particular, the decay rate of Z — > /it in the 
SUSYE331 model can be written as follows: 



9z m z 
24tt 



4mn V2 



m| 



[a] + b] 



m 2 z / 



+ 6a E b E — 



rrir 



(28) 
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where 



- 2 - '-^ and /„ 0,4- 4- 



0,£ = Tf — QtS^r — 



For the decay rate of Z — > /ir, because the form of effective operators of Z\xt 
in the considered model is the same as that of the MSSM, the form of decay 
rate of Z — >■ /it is the same as that established in [9]. Especially the result is 
given as follows: 



where 



9z m z 
2An 



+ 



R 



+ 



rn 



D1 



m z 



2 1 



rn 



m z 



U R 



,(29) 



L,R 



\Z (~iZ 
L,R L,R- 



Comparing the two results taken from Eqs.(28)and (29), we obtain the re- 
lationship between the two branching ratios corresponding to two processes 
Z — >■ 11 and Z — > fir. Our calculation is consistent with results given in [9] 
such as: 



BR(Z ->■ fir) 



= crrir 



+ 



R 



2 1 

+ 2 



m z 



2 1 
+ - 



D 



m z 



R 



BR(Z ->£+£"), (30) 



where c = (a 2 + bj)~ l = (l/4-s 2 v +2s$ v )- 1 ~ 7.9 and BR(Z ->■ £+£") ~ 3.4%. 



5.^.5 Z' ->■ /xr 

Similar to the case of the Z boson, the decay rate of Z' — > £ + £~ can be 
determined in the below formula: 



2\ 1/2 



24vr 



(a? + bf) ( 1 - ^ ] + Qa' e b' 



9z' m Z' ( 2 , r 2\ 

where af + 6f = aj + 6^ = 1 /c and 

4 = (4c 2 w - l) (T| - Q e ) + 3c^X, = -a 

bf e = - (Ac 2 w - l) Q £fl + 3c 2 w X eR = -b e . 



ml, 



1 e m 2 
m z , 



(31) 
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The decay rate Z' — > fir is similar to that of Z — > fir. It leads to the relation 
between two branchings as follows: 



BR(Z' fir) = c m%, 



+ 



+ 2 



xBR(Z'^ft), 



-D 



z' 



m z > 



1 

+ 2 



rn 



m z 



-D z ' 

U R 



(32) 



where 



,Z'(2) 



L,R - _,2 L,R -T L,R 

Illy, 



n z ' 

°L,R- 



(33) 



3.4-4 r — > f 1 ^^ 

In the SUSYE331, the effective Lagrangian described the r — > fififi decay 
can be deduced from the effective operators given in Eq.(26) combining with 
those induced by the effective operators of lit with Z or Z' or photon as well 
as Higgs boson. The general study was presented in [9]. The contributions 
from the box diagrams and vector boson exchanges to the branching of the 
considered decay rate are sub-leading ones even when tan 7 is large or small. 
However the contribution from the Higgs bosons exchange to decay rate is 
large for large tan 7 (the interested reader can see in [20]). Hence, in this work 
we will split each type of contributions to the considered branching ratio. First, 
let us consider the case of absence of Higgs exchange: the effective Lagrangian 
can be deduced from the effective operators given in Eqs.(15), (18), (19), (25) 
and (26). The explicit formula of effective Lagrangian is 



-.cff 



+ (/i W) (Fjfp*^ + i^V<V*0] 
+ 2e 2 (Dl^^ + Dl^r)^ 
x (pa^fx + /iV M /I c ) + h.c, 



(34) 



where 



>Z> _ m Z A Z'(1) A Z'(2) 

l L(R) ~ m 1 L,R "I" L,R , 
1 •I' 71 



(35) 



F L(R) - B L{R) + 29z c 2W^L(R) ~ ^Qz'^W ^L(R) ~ ^^{R)-, 



^L(R) — Bl(R) + 9z s W^L(R) 9z' s W^L(R) e ' 2( ^L(R)- 



(36) 
(37) 
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Here we also assume that, as in the case of MSSM, we ignore the contributions 
°f ^lri Clm ^lr an d D% R to the above effective Lagrangian. This leads 
to the branching ratios of decay r — > fififi is 



BR(r — > ji ji + ji 



8G F 



2\F£ L \ 2 + \F£ R \ 2 +\F^ L \ 2 + 2\F^ R \ 2 



+ 4e 2 Re (dI (2Fj? + F?) + Dl (FJ? + 2Fj*)) 
+ ^(|^ + |^)(log^-^ 

x BR(r _ ->■ fTv^Vr). (38) 



3.5 HfiT contribution to r — >■ \x\x\x 



Contribution of Higgs exchange in the SUSYE331 model was investigated in 
[20], where the corresponding effective Lagrangian for this, is given by: 



x(A£/iT c + A£//T) + H.c., (39) 



where 



Cst, (^- + ^-)- (40) 

Noting that factor (m^tan7) can not be ignored if value of tan 7 is large 
enough. This factor causes a shift to F% R and F^ L , explicitly 

8F£ R = ^2G F m^m T CA R , 5F£ L = y/2G F m^m T C^L- (41) 

The individual contribution from Higgs exchange to BR(r — > fi"fi + fi~) now 
is: 



Km r C) 2 (|A£| 2 ) + |A£j 



2 



BR(r -> fi fj+fi = ^ — BR(r ->■ fi ^ T ).(42) 



This contribution will be suppressed in the case of small tan 7. 
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4 Numerical results 



In this section, let us numerically study the cLFV decays of the tau lepton 
r — > p'j, t — > ppp and Z — > pr. For this purpose, we first study some 
constraints on values of parameter space in the SUSYE331 model. 



4-1 Implication on the parameter space in SUSYE331 model 

In this section, we pay attention to discussing on constraint of parameter space 
caused by experimental cLFV bounds (l)-(3). As mentioned, this topic has 
been carefully studied in many models beyond SM. Especially ref.[9] not only 
indicated regions of parameter space satisfying experimental bounds but also 
discussed in details the correlation between dipole and non-dipole kinds of 
contributions to cLFV decays in each region. The most important assumption 
here is that sources of cLFV come from only the mass terms of sleptons in the 
soft term, namely only left-and right-handed slepton mass matrices have large 
p — t entries. Let us compare the sources of cLFV appearing in the MSSM 
with those of the SUSYE331. Because of the absence of the right-handed neu- 
trinos, the MSSM contains only three sources of cLFV, particularly LFV in 
left-handed charged sleptons, left-handed sneutrinos and right-handed charged 
leptons. But the left-handed sleptons and their sneurinos live in the same dou- 
blet of SU(2) L in the MSSM, the origin of cLFV in the left-handed charged 
sleptons and left-handed sneutrinos sectors are the same As a consequence, in 
the MSSM there are only two independent sources of cLFV. Due to the ap- 
pearance of the right-handed neutrinos in the SUSYE331, there also appears 
one more source of LFV. Furthermore, there exist two Higgs multiplets in the 
model, i.e triplet p and antitriplet p' which independently generate masses of 
charged and neutral sleptons at tree level. These two Higgs multiplets also cre- 
ate at least two new corresponding mass terms of sleptons in the soft term, as 
shown in details in [19,24,20]. In general, there are at least four independent 
sources causing the cLFV in the SUSYE331 where each source is parame- 
terized by a mixing angle defined in the last subsection of the appendix A. 
This parameterization creates the similarity between the SUSYE331 and the 
MSSM. We exploit this advantage to make some prediction for the SUSYE331 
basing previous investigation of the MSSM. For example, if all of these sources 
are appeared, the branching ratio of decay r — > p'-f will be much enhanced 
than that in the MSSM, even it could greatly exceed the bound of experiment 
given in (1). As considered in the MSSM, the existence of maximal LFV mix- 
ing in the left-handed slepton sector means the LFV mixing in both charged 
and neutral sleptons, and their contributions to cLFV decays are much larger 
than contributions caused by right-handed LFV mixing. In addition, if the 
superparticle spectrum is relative light, namely the parameter space are set 
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below 1 TeV, there will exist very small regions of parameter space satisfying 
experimental results. The similar situation also happens for the SUSYE331. 
However, in the SUSYE331 four cLFV sources are independent then there 
exists a situation that the model contains only one left-handed maximal LFV 
source while others vanish. It is easy to realize that in this case, the values of 
branching ratios of cLFV processes in the SUSYE331 are smaller than those 
of the MSSM and the satisfied regions of parameter space will be wider in the 
scale of (9(100) GeV. On the other hand, if four LFV sources are presented, 
the predicted results of the considered model consistent with the experimen- 
tal results at the TeV scale. In the next subsections, we will examine the 
influence of LFV sources appearing in the soft term of the SUSYE331 on the 
cLFV decays of the tau and the Z boson. In the numerical investigation, we 
just pay attention to the case of soft parameters at 0(100) GeV scale because 
this scale allows the existence of small values of slepton masses which can be 
detected by present colliders. In fact, the detailed investigation to determine 
different properties of cLFV branching ratios among different regions of pa- 
rameter space is really needed, as done in many known SUSY models. In the 
SUSYE331, this work is more complicated because of the addition of many 
new particles so we will come back this interesting topic in another work. 

We would like to emphasize that the sleptons gain mass through main sources 
which come from soft terms and interacting terms of sleptons with Higgs 
through F-and D-terms. As mentioned in [25], the D-term gives contribution 
to mass of slepton that contribution is proportional to the quantity {w 2 — w' 2 ) 
while F-term, gives contribution to mass of slepton that is proportional to 
(\ a w) or (\ a w'). Because the VEVs u and cu' break down SU(3) L to SU(2) L 
so these values can be set to the TeV scale. However, the SUSYE331 contains 
the Tachyon fields, the removal Tachyon fields leads to a condition [24]: 

\w 2 -w ,2 \~ \v 2 -v' 2 \< 246 2 (GeV 2 ). (43) 



The conditions given in (43) show that the D-term generates mass for slepton 
at the 0(100) GeV. Furthermore, if the R-parity is assumed to be preserved, 
the coupling A a must be zero. It means that the contributions from both F- and 
D-terms to masses of sleptons all in the range of hundreds GeV. So if slepton 
masses are in the range of 0(1) TeV scale, the dominant contributions are 
from the soft term, as in the MSSM. 

Next, let us discuss on \x p and fi x which play a same role as fi parameter in 
the MSSM. According to [19,24], the requirement of the canceling all linear 
terms at tree level in the Higgs potential leads to conditions: 



x 3 - 4s 2 



w 



2 ( 12 2 i 12 2 \ i i 2 I2\ 

c w [W - W ' + U - U j + - \ V 1 - V J 
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4 + H = - 2 (3-4 s * w ) K W ' 2 ~ w2 + U ' 2 ~ u2 ) +2 ( v2 - V ' 2 )] ■ (44) 
They result in two consequences: 



K + ±<\< 2 f^ ] M 2 w and |^ + 4m;|<3-^-J<, (45) 
where M w ~ 80 GeV is the mass of the W boson. 

If all parameters fi 2 ,, fx 2 , m 2 x and m 2 are positive, the constrains given in (45) 
set values of them in the scale of (9(100) GeV. 

It is known that three branching ratios r — > /ry, r — > /i^u/i and Z — > /it 
depend on some of following quantities: A z ™ , B%$ , Cf R f\nd Df/'\ 
In the different regions of parameter space, where some of these quantities give 
dominant contribution, there are precise correlations among branching ratios. 
In particular, with three considering branching ratios, we have two cases as 
listed below: 

D 1 - dominance: It is easy to get the relation concerned in [9] 

BR( ; ^4 - 2.2 x ID" 3 , (46) 
BR(r^/i 7 ) ' V ; 

and a consequence is 

BR(r -> fififi) < 10- 10 . (47) 
A z - dominance: In this case, we get 

/, 



BR(r — > umi) 



0.53 
0.67 



BR(Z ->■ fir) (48) 



Note that in the SUSYE331 model, although ^l(r) receives contribution 
from diagrams which are similar to those of Af^, A R ~ because the 
ty — u)/u>' ~ 1 (c2 7 ' = 0) and Af depends on only diagrams containing 
right-handed sneutrinos. 



In next subsection, we consider numerical computation for cLFV decays of 
the tauon and the Z boson. Remember that the effective couplings used for 
numerical studying are established in appendices B, C, D and E. The standard 
loop integrals are given in [9,23]. 
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4-2 In the case of small tan 7 and light mass spectrum 



4.2.1 t — y fj,^ 

Let us consider the numerical studying of branching of r — > /ry decay. The 
analytical result given in Eq.(27) depends on the effective coupling D 1 ^^ 

which can be divided into three parts, -D^r) = + Dlffy + D]^. The 

analytical formulas of Dj^L, ^Iur) anc ^ are §i ven i n appendix B in which 

only -D2(R) does not depend on tan 7. In addition, from the experimental 
bound (1), we can obtain the constraint on the effective couplings, namely 
|£>Z >H | < 2.5 x lCr 9 [GeV 2 ]. 

We would like to remind that the diagrams which contribute to the D\ are 
collected from three LFV sources: left-handed charged slepton, left-handed 
sneutrinos and right-handed sneutrinos sectors, while the diagrams contribut- 
ing to D\ are only collected from the charged right-handed slepton sector. 
So the values of D\ are predicted much larger than those of D\. Another 
thing we want to emphasize that since the SUSYE331 has many additional 
particles, we expected that the additional particles can modified the predicted 
results of cLFV decay in the model under consideration. For more details, let 
us consider numerical studying in this decay mode. 

First we study the effects of LFV sources on D\ in the case of small tan 7 
as well as the presence of all of the three left-handed LFV sources, especially 
we fix tan 7 = 3 and 6 L = 9 H = 6 Pr = tt/4. As in the MSSM, Df b) is 
dominant contribution to D\. Fig.l shows the dependence of D]}® on soft 
parameters m~ Lz = = and m x , while other parameters are fixed. 

The predicted results given in left panel of fig.l fully consistent with the 
experimental results if the domain of parameter is close to the value of 
mi 2 and all soft parameters are set at (9(100) GeV. We also remind that 
mixing mass term between left and right slepton, , is small if the model 
under consideration has maximal mixing sources. However if ji p is set at TeV, 

the domain of parameter m^ 2 , where the values of D 1 ^ match experimental 
bound, is more extensive. For more details, the reader can see in the right 
panel of fig.l. 

Let us consider the case m~ LZ = 1 TeV. The results given in fig. 2 predict 
that the values of D]^ exceed the experimental bound if the remaining pa- 
rameters are assumed in region of (9(100) GeV. Even if the values of fi p are 
assumed as large as 1 TeV, the predicted value of is larger than the 

experimental bound. In this region of parameter space, the predicted results 
in the SUSYE331 are much similar to those in the MSSM [9]. As mentioned 
in [9], the large values of A T as well as small values of ji p and masses of slep- 
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Fig. 1. isocontours with tan 7 = 3., = m u L3 = m i>R 3 an d 

m^ 2 = mi> L2 = m<u R2 = 300 GeV, 9l = 0u L = &u R = and fi p = 140 GeV 
(ITeV) in the left (right) panel. The solid and dashed lines correspond to = 300 
GeV andm B = -300 GeV. 

tons must be required to keep value of the total |-D 7 | below the experimental 
bound. However, the fig. 2 displays the predicted values of D\ are much greater 
than those in the MSSM. It means that in order to obtain the values of D 1 ^ 
being consistent with experimental bound, the values of A T predicted in the 
SUSYE331 should be much larger than those in the MSSM If in the limit of 
large values of A T , the model leads to the appearance of Tachyon sleptons. 
It means that the model under consideration does not contain the region of 
parameter space such that there exists a large difference between the values 
of slepton masses. The fig. 2 also shows that the values of D\ exceed to the 
experimental results when parameters expands into range of TeV. 



D£ IW [10~* GeV -2 ] Dl""[l0"' > GeV -2 ] 




100 150 200 250 300 100 150 200 250 300 

mi[GeV] mj[GeV] 



Fig. 2. L>2 isocontours with tan 7 = 3., = mp L = fnv R2 and 

m^ 2 = Tay L2 = rriv R2 = 1 TeV, 9l = 0u L = @u R = vr/4 and fi p = 140. GeV 
(ITeV) in the left (right) panel. The solid and dashed lines correspond to ms = 300 
GeV and m B = -300 GeV. 

So in the SUSYE331 model with the case of the maximal mixing happening in 
all three sources (left-slepton, left-sneutrino and right-sneutrino sector), the 
scale of all slepton masses should be the same order, in range of TeV or in 
range of (100) GeV. 
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Now we consider another situation that happens only in the SUSYE331, not 
in the MSSM. It is the case of only one left-handed LFV source appearing in 
the model. Looking at analytical formulas of effective couplings, three sources 
contributing to left-handed effective coupling are parameterized by three mix- 
ing angels 9l, Qv l and 9& R . Two of them, 9a L and 9c R , relate with diagrams 
containing sneutrinos propagators while the remain relates to diagrams with 
charged sleptons propagators. Numerical computation indicates that if the 
mass spectrum of superpartner particles are in the range of (9(100) GeV, con- 
tribution from sneutrinos exchanges are larger than those of charged slepton 
ones. So if two sneutrinos mixing angles vanish, there only one source of left- 
handed mixing 9l which generates relatively small effective couplings. The 
experimental bounds then are easily satisfied even in regions of light mass 
spectrum. Our numerical investigation will focus on this case. In particular, 
mixing angle parameters are fixed as 9 R = 9„ L = 9„ R = and 6 L = ir/4. 
The Fig. 3 displays as function of the A T and fi p while others are fixed: 
mi 2 = 1 TeV, m R2 = = mp L2 = rni,^ = m u R2 = m v Rz = m R - The results 
given in the fig. 3 illustrates that in the considered limit, we can find the region 
of the small absolute values of A T in which we can obtain the values of \Dl\ 
satisfying the experimental bound, particularly: 

- The value of tub should be smaller than that of m\. 

- If the value of is closer to that of m\, the parameter space of A T and, \x p 
satisfying the experimental bound of \D1\ has been expended. 

The predicted results given in fig. 3 show that if in the SUSYE331, only one 
source of lepton number violation in the charged slepton sector, the model can 
contain the region of parameter space in which the all soft parameters are set 
to 0(100) GeV except m^ 2 is set to TeV. In this region of parameter space, 
the predicted results on the r — >■ /ry are matched the experimental bound on 
the decay. The existence of the soft parameter space below TeV scale leads to 
hope that the sleptons can be detected by LHC. 

Finally, we consider the LFV effect in the right-handed charged slepton sector 
to the LFV in the tauon decay. We assume that there is only the maximal 
mixing of right-handed slepton, i.e., s R = c R = l/y/2 and all other mixing 
sources are set to be equal to zero. Under this assumption the Feynman di- 
agrams in the model under consideration contributing to D 1 are exactly the 
same as those in the MSSM [9]. In the fig.4, absolute values of D R {a) and D R {b) 
are rather small, even a bit smaller than those in the MSSM. From this inves- 
tigation, we see that A T should have the same order with the mass parameters 
of sleptons. 
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Fig. 3. isocontours with tan 7 = 3., m i 2 = 1 TeV, 

6 L = vr/4, Or = 6} L = 9 Vr = 0, = (only LFV in {rh L ,f L } sector). 

For illustrations, we choose three choices of parameter space (mg, m\, rn^^, m R ) 
in GeV: (200, 250, 300, 200) ( solid) , ( 100, 300, 100, 200) ( dashed), ( 100, 300, 
300, 100) (dotted). For each example, the center line corresponds to the value of 
Z?2 = 0) two other ones limit the region where \Dj^\ < 2.5 x 10~ 9 [GeV -2 ]. 
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Fig. 4. Isocontours of D R ^ (left panel) and isocontours of Dj^ (right panel). The 
parameters are tan7 = 3., m^ 2 = 1 TeV, Ol = @u L = 0p R = 0, Or = 7r/4 and 
Hp = 150 GeV . 

4.2.2 Correlations 

Next, we consider the process r — > /i/i/i. The amplitude consists of the contri- 
butions from the effective couplings Al(r), B l ^,Cl(r) and -Dl(r) which are 
denoted the coupling of Z, Z' gauge bosons, photon and Higgs to charged lep- 
tons Hl(r), t l(r). Each type of diagrams picking up the each particle exchange 
contributing to the LFV in the considered process depends on the region of 
parameter space. In next work, we will study the contribution of each effective 
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coupling to the LFV of the r — > /i/i/z process in two typical cases. 

Only maximal mixing in the charged slepton (jl, f) 

Let us consider the process r — > fififi in the limit of small tan 7 and in the 
case of existing only maximal mixing in the charged slepton, sl = cl = 

and sji = Sy L = s„ R = 0. This constraint leads to Ajf' = Aj^' = 0. As 
mentioned in the last part, it is interesting to consider the region of param- 
eter at 0(100) [GeV]. One can check that Af^ is dominated by Af^ and 
total effective couplings contributing to the branching ratio of the considered 
process are -DLm, ^l(R)' an d B^ L(R) . Indeed, looking at the analytical 

expression of these couplings we found that only the analytical expression of 
D\ is affected by left-handed slepton parameters such as: A T , A^ T . Hence, in 
order to look for the effective coupling giving dominant contribution to the 
considered branching ratio, we will study numerically two cases: 

- The A T , A^r parameters are fixed and the other soft parameters are changed. 

- The soft parameters are fixed and the A T , A^ T parameters are changed. 

First, we assume that A T = A^ T = 0. To assess the contribution of the each 
effective coupling Af and D\ into the BR (r _ — >■ fi~fi + fi~), we define two 
factors f^z and fr>i such as: 



9% 




(l C 2W + S w] 


+ 


A Z R 


(j c 2W + ^ s w) 




2\F£ L \ 2 + |FH 2 


+ 1^ 


2 + 2\F£ R \ 2 



(49) 



and 



8e 4 



fm = 



+ 



2 ) h (i) - 



MS 



(50) 



where 



MS = 8e 4 (\Dl\ 2 + \Dl\ 2 ) 



+2 \Fr L \ 2 + \Ft r \ 2 + IF; 




M|2 



(51) 



The factors /^z and f^,-, given in Eqs. (49), (50) quantitatively measure con- 
tributions of ^.f(fl) and -D^m to the factor |F ML ( K )| and the total branching of 
(t~ — > [i~ ji + fi~) decay, respectively. 
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Fig. 5. Correlations among Af, F^ L< ~ R) and D[ with A T = 0. The contours f A z, f D i 
and BR(t — > ^7) are denoted by dashed, dotted and solid black lines. For illus- 
trations, numerical values for parameter space (tub, rn\,'rni 2 ) are chosen as 
(100, 100, 1000, 100) [GeV] (left-panel) and (100,300,1000,100) [GeV] (right panel). 

Looking at the Eqs.(37), (36), the coefficient F^^ } depends on the other 
factors. If A z gives dominant contribution to F^^ ] , it is convenient to define 
regions where the factor f^z satisfies 1.05 > f^z > 0.95 as the A z -domination 
regions (over F^^). On the other hand, we also make convention that if 
fm < 0.05 then the dominant contribution to the branching ratio of r — > 3fi 
is given by F— type of couplings and if 1.05 > foi > 0.95 then the dominant 
contribution to the branching ratio of r — > 3// is given by D— type of couplings. 

The current experimental upper bound on the branching ratio of the process 
BR(t — > /ry) is smaller than 4.4 x 10~ 8 . The results of our calculation to this 
process given in fig. 5. The experimental bound of the branching is denoted by 
the solid black lines. One can see that this process satisfies the experimental 
bound in regions of large parameter space of /i and m L ~ 3 . 

The regions of parameter space where Af gives dominant contribution to 
branching ratio r — > 3/i must satisfy both conditions: 1.05 > \/a z \ > 0.95 and 
I/D7 1 < 0.05. One can check in fig. 5, there is very narrow region of /i p and 
m^ 2 parameter space satisfying both conditions, particularly the parameter 
li p has very strong constraint \i p -r [250 — 280] GeV. The fig. 5 also shows that 
if D 1 gives dominant contribution to the considered decay mode, /i p fluctuate 
slightly around the value of 300 GeV. 

Let us consider the effects of A T on the branching of the considered decay 
mode. We would like to emphasize that A z does not depend on the A T while 
the value of l-D^I depends on the sign as well as amplitude of A T . In particularly 
D]} ^ is proportional to (A T + |/i p tan7), the values of A T will affect on soft 
parameter space region where D\ gives dominant contribution to the BR(r — > 
3/i). These regions of soft parameter can be larger than that in the case of 
A T = 0. 
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Fig. 6. Branching ratios Z — )• [it (left panel) and r — > 3fi (right panel) as functions 
of mg, Three numerical values for parameter space (m\, fi p , wi^ 2 , m i 3 i 771 _r) [GeV] 
are chosen: (100, 100, 1000, 100, 100)-black line, (200, 100, 1000, 100, 100)-green line, 
(100, 300, 1000, 100, 100)- blue line. The red line in the right panel present the bound 
of experiment BR(r — > 3fi) = 2.1 x 10 -8 . 

Now let us consider Z — > fir and r — > fififi decays. The branching ratios 
of decays Z — > \ir and t~ — > fififi are presented in (29) and (48). The pre- 
dicted branching is shown in fig.6. The numerical branching ratio for Z — >■ fir 
has the maximum of 10 -9 for the soft parameters setting to (9(100) GeV. 
This predicted result is very suppressed with the present experimental bound. 
However, in the same region of soft parameter space, the predicted result for 
t — > fififi can reach to the experimental result and even one can find some 
regions of parameter space that predicted our result exceed the experimen- 
tal bound. One can check in the left panel of the fig.6, the branching values 
for Z — » fir hardly change when we change the value of the parameter mg. 
Moreover, as we increase the values of the soft slepton mass parameters these 
branching values decreased. Therefore, to increase the value of the branching 
of Z — >• /ir, we have to change the parameter \i p . 

The fig. 7 shows the values of branching ratios of the decays in the plane 
tub — fi p . In this case, we choose A T = and other parameters are chosen 
so that the experimental branching decay of BR(r — > fififi) is satisfied. We 
can see that the bounded regions of BR(r — > fij) supports the small values of 
both remain decays. In the case A T ^ 0, because only BR(r — > n^) depends 
on the values of A T and as we have shown in the above section, there will 
exist a possibility where D\ vanishes. This case allows BR(t~ — > fi~fi + fi~) 
can reach the limits of experiment of order O(10~ 8 ) whilst BR(Z — > /xr) is 
still in maximal order of O(l0~ 9 ). 

Only maximal mixing in the right-handed charged slepton sector 

/i, T. 

Now we come to consider another case, only maximal LFV in right-handed 
sector of charged sleptons (sr = l/y/2, Sl = Sv L = Sy R = 0) where re- 
gions of parameter space can be available in range of 0(100) [GeV]. Both the 
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Fig. 7. Contour plots for branching ratios r — > /j, /i + ji (dot- 
ted), Z — > fir (dashed) and r — > yj (solid) with A,- = and 
(m A , m Za , m £s , m^) = (100, 150, 1000, 100, 300). 

SUSYE331 and the MSSM models are similar to each other in this case. So 
we just discuss more on correlation among effective couplings. As shown in 
the left panel of fig. 8, the bound of experiment of BR(r — > y/y) rules out the 
large values of BR(r — > 3//). This leads to the result BR(r — > 3fi) < 10~ 9 if 
A T = 0. The right panel shows that in the case of only large LFV in right- 
handed charged slepton sector, the BRij — > 3y) is in maximal order of 10 -9 , 
even in the case of non-vanishing A T and A T makes Djf suppressed. For the 
BR(Z — > fir), this case is much smaller than the previous case. Also we can 
see a difference from the case of pure large LFV in left-handed charged slepton 
sector: the D 7 -domination regions now lies on the small values of y p while the 
large values of /i p are ruled out by the condition BR(r — > fij) < 4.4 x 10 -8 . 



5 Conclusions 



In present paper, we have studied the LFV decays of the tauon and the Z 
boson in framework of the SUSYE331 model, and have mainly focused on 
two-generations slepton mixing, namely both lef- and right-handed y — r slep- 
ton mixings. In order to obtain the relevant diagrams, we have combined the 
mixing of sleptons, that of charginos, Higgsinos, gauginos as well as the inter- 
actions of gauge bosons with leptons and the Yukawa interactions. From these 
diagrams, we obtained the effective operators relating to the considered cLFV 
decays. This leads to the analytical expression of the branching ratio of the 
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Fig. 8. Contours in [i p — rn^ plan (left panel) and plots of branching ratio of 
r — t- 3/i (right panel) in the case of tan 7 = 3 and A T = 0. The respective con- 
tours are BR(r — > fij) (solid lines), fjj-y (dotted lines) and BR(r — > 3//) (dashed 
lines) with numerical values of parameters {mB,rn^,mj^ 2 ) = (100,100,1000) 
(m^ 2 = m^ 3 = m^). For the plot of BR(r — > 3/i) four choices of parameter space 
(tub, /x p ,m^ 2 ,m^ 3 ) are: (100,100,1000,100) (black), (200,100,1000,100) (green), 
(100, 200, 1000, 100 3 ) (blue) and (100,300,1000,100) (red). 

considered decay processes and the contacted relation between decay rate of 
non LFV decay mode with that of cLFV decay mode. Our analysis is carried 
out in the limit of small tan 7. The detailed predictions in our model depend 
strongly on the SUSY parameters and left- and right-handed slepton mixing. 
Consequently, we have firstly considered the effects of SUSY parameters and 
the mixing of left- and right-handed sleptons on the r — > fij decay such as 

• In the case of the maximal LFV mixing, the mixing mass terms between 
left-and right-handed sleptons {m% t , m R t ) are small, our results are only 
consistent with the experimental bounds if the domains of parameter m^ 3 
are close to those of m^ 2 whenever they are set to the TeV or 0(100) GeV 
scales. It means that in the case of maximal LFV mixing, the slepton mass 
parameters are in the same order. 

• If there is only the LFV in the charged left-handed sleptons sector, we can 
find some regions of parameter space that allow above cLFV branching 
ratios matching with the experimental bounds. Especially the slepton mass 
mi 3 is set at 1 TeV while the other, m^ 2 , is set at (9(100) GeV. Noting 
that the value of parameter tub should be close to that of m\ and if the 
value of is closer to that of m\, the parameter space of A T has been 
more expanded. In the case of left-handed LFV sector, one important result 
deduced from our numerical investigation is that although the SUSYE331 
model contains much more supersymmetric particles as well as LFV sources 
than the MSSM, there still exist of some wide regions of parameter space 
which allow not only masses of sleptons but also fi p keep the small values 
enough to be detected by colliders. 

• If there is only the LFV in the charged right-handed slepton sector, in order 
to match the experimental bound, the value of A r should be the same order 
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as that of other soft parameters. Our result is similar to the predicted result 
in the MSSM. 

Based on the parameter space satisfying the experimental bound on r — > /ry 
decay rate, we consider the branching ratios of r — > fififi and Z — > fir. We have 
concentrated on the LFV in the charged left-handed slepton sector with only 
8l = 7r/4. In this case, the dominant contribution to the considered decay 
mode comes from Az if fi parameter has very narrow constraints such as 
fi E (250 4- 280) GeV, while D 1 gives dominant contribution to the considered 
decay mode if the \i parameter has a small fluctuation around 300 GeV. If we 
set the value of A T = 0, the constraint on the fi parameter can be expended. 
Similarly, by numerical study on the branching ratio of r — > fifi/i decay mode 
in the case there exists LFV only in the charged right-handed slepton sector. 
The small value of fi p giving the dominant contribution to considered decay 
mode coming from D 1 was obtained. In both cases, our predicted results of 
Z — > ji'j are very suppressed. 
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A Interacting Lagrangian and notations 



Some of the interacting vertices used in this work was given in [20]. In this 
section we list the rest part of interacting Lagrangian which is necessary for 
completion our calculation. Here we use the relation = (£ C L )* where £ C L is 
the superpartner of a lepton £ C L in the model. 



A.l Interaction between left-right slepton sector with neutral higgs bosson 



These terms come from two sources: 
(1) From F-terms: 



J L p p° (Y^lf* + Y T r* L rr) + H.c. 
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(2) and from soft breaking term: 

^UhO = ~habLaLP IbL + H.C. 

= -K b (v aL p'{ + l aL p'° + V c aL p' 2 ) l c bL + H.c. 
-Y,A^°~p L ~p c L - Y T A T p'°f L T c L 

-Y T A L ^p L rl - Y T A* T p'°f L jli + H.c., (A.l) 
where we use new notations that is identified with those in [9]: 

^22 = h-'ii/j, = Y^A^, ^33 = h TT = Y T A T , 
= = Y T A^ T , h 32 = h Tfl = Y T A^ T . 

The total interaction part of (H°l c L l L ) interactions is: 

4£o = -Y, (i/i P p * + A^ p L p c L - Y T (±p pP °* + A T p'^ r L ft 

-Y T A^p'°p L r L - Y T A* T p' f L pi + H.c. (A.2) 

A. 2 Gauge boson interactions 

This kind of vertex is only contained in gauge invariant kinetics of all fields in 
the theory. In this work, we just study on cLFV in lepton sector so the related 
part of the Lagrangian is [19]: 

Akinetic = (D^DuP + (L>V) f ^y + fp^D^p + ipcPDj 

+ (D^LiL^D^La + iLiL^D^LiL 
+ (D 1 ,l c iL yD 1 ^ iL + il^D^ 

- \F^F a ^ v - \f^F^ v + iX^Dfa + f\ B a^X B (A.3) 

where % = 1,2,3 is family index, a = 1,2, ..,8 corresponds to eight gauge 
bosons of SU(3)l group. Covariant derivatives D^, and correspond 

to triplets, anti-triplets SU(3) L , singlet SU(3) L and adjoint presentation of 
SU(3) L . They are defined as follows: 

F a{ iv = dfjVav — d v V a n — gf abc V blj y ci/ , F^ = d^By — d v B^ 
D^ = d^ + igT a V ail + tg'XT 9 B^, 
D^ = d,-igT a *V a , + ig'XT % B^ 
D 1 ^ = d, + ig'XT 9 B„ 
Df;\ a v = d II \ a v -gr bc V b \ c v . (A.4) 
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Here X denotes U(l) hypercharge, f abc is structure constant of SU(3), T 9 is 
the generator of U(l)x which is defined by T 9 = 1/V6 diagonal(l, 1,1). We 
just pay attention to neutral bosons in covariant derivatives so they can be 
written as [19,16]. 



= d» + ig (T 3 IV + T 8 V SlI + tT^XB^ + T% + T 5 V 5fl ) , 
D N = d -iV nc 

= 8,-ig (T 3 VV + T 8 V 8 , - tT^XB,) - ig (r% - T 5 V 5/l ) , 
D^d^ + ig'XT^B^ 

d l ^v = d,\ a v - g (r 3 % 3 + r% 8 + r A x + f a5 X) X °v ( A -5) 

where gauge bosons W 3 , W 8 , W4, and B relate with physical states according 
to the transformation: 



fw 3 \ 

B 
\W 4/ 



s w 

^3 



C^Cg'Cw 

C ip K3 — S v KlK2 

tyy( c <p K l+ s ip K 2) 

\fZCnt 



\JHcwCqI 
twi s <p K l — c <p K 2) 



\ -tg,(c ip K 2 - S^Ki) -tg,(s^K 2 + C^Ki) 



Sd'Cw 


«2 



'a" 
z 

Z' 



(A.6) 



where some new notations are used: 



to = tan# = — , t 2 e = tan(20), sq. 



w " c w ^l + &l e 

3\^2sw nr? 7 3\/2sw 



g 6^is w _ r— 2 - 



? 2 



K 2 = \Jl- 4:S 2 6 ,c w , k 3 = s w - Sc w sj,. (A. 7) 

In the SUSYE331 model, we have all 9, <p and 6'<Cl. Thus, we can use the 
approximation sin# = sirup = sinO' = tan# = tamp = tan#' = to simplify 
the calculation. We also take the approximation: 

K 2 ~ 1, K 3 ~ S w . 
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In addition, W$ and W' A make of a physical neutral non- Hermit ion gauge boson 
X° which is defined by the combination: 



(A. 



So we can rewrite the above covariance derivatives in the form below: 



Dp —d^ + xeQAp + ig z (t 3 - s 2 w Q) 
+ igz> [(4c 2 ^ - 1)(T 3 -Q) + 3c 2 w x] . 
f ieQAp + ig z (-T 3 - s 2 w Q) Z i 



+ igz> 



(Ac 2 w -l)(-T 3 -Q) + 3c 2 w x]z^ 



D in -9^ + ieQAp - ig z s w QZ^ 
+ ig z , [-{Ac 2 w - 1)Q + 3c 2 w X 

where we have defined 

gc<p 



z'„ 



gz 



CwCgi 



— and 



gz> 



9 



gc^ 



(A.9) 



For the charged gauginos, we have: 



± 



W* = -^-^, F^^p. (A.10) 



This leads to the covariant derivative of charged gauginos: 

D^W* ~ dpW* ± % (eAp + gcwZp) W ± = d^W* ± % (eA^ + g z c 2 w Zp) W 

= dpW* + iQ w (eAp + g z c 2 w Zp) W ± , 
D L Y± = d,Y± ± i (eA, + g ^w + 2^4) + s v K lK2 

S v (C 2 W + ZCwS 2 ,,) - C^KiK2 z , \ y± 

9 2c w c e > 7 

-d^ + iQy* (eA li + ^g z c 2W Z li -^g zll 4Z' l ^Y ± . (A.ll) 

From these two formulas, we can deduce the vertices of neutral gauge boson- 
charged gaugino-charged gaugino . 
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A. 3 Gauge boson-slepton-slepton interactions 



This kind of vertex comes from the part [19]: 



hiv ~ 2 



+ 



\/6 



I r 



ya 



B„ 



where i = 1,2,3 is the flavor index and a = 1,2, ..8 is generator index of 
SU(3). For the {/2,t} sector with neutral boson we have: 



'-uv — 



tg 



W 
+ 2 

ig 

2 



c w 
-2-SwA^ 
2c w 



K\Cw 

cw 



x [d il v T v T — i> T d^i> T ) 



ry . C 2 W ry, 



KlC W 



-Z' 



+ 



e tW ry, 



i[eA^-e twZ^ + (>rV - r c ^f c 



— z 



— i 



+ 



-g z Z^ + -gz'CawZp 

eAf, + ^g z c 2 wZ fl - ^g z >c 2 wZ' fl 

'gz'C w Z^] (d»v\v\ - z> r cKV\) + (r ->■ y) 
(eA, - gzs^Z, + gz>s w Z'^ (>rV - r^r 



+ (r° -f J* c )] 



(A.12) 



p 



V 



V 



Fig. A.l. Notations of directions of scalars and fermions. Here denotes a photon 
A, Z or Z' boson. 



Interaction vertices of photon, Z and Z' bosons relating with our calculation 
are summarized in tables A.l, A. 2 and A. 3, respectively. We denote directions 
of momentums in fig.A.l. For simplicity, we omit spinor index in the formulas 
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of boson-fermion-fermion vertices. The precise formulas of this kind of vertices 
is easily deduced using rules concerned in [11]. 



Table A.l 
Photon vertices. 



vertex 


factor 


vertex 


factor 


Photon-scalar-scalar 
(scalar ip : H, f ) 


— ieO.Jv + v'Y 






Photon-spinor-spinor 
(spinor tp: fermion, Higgsino) 








Photon-boson-boson 

W +P W~^A U 


ie\p +p , p-n, pau] 


Y+py~v-A u 


-ie\p +p , pau] 


Photon-Higgs-gauge boson 
A^W u pi 




A^Y u p 2 


\{ieg)g^v 


Photon-gaugino-gaugino 
W+A^W + 
Y+A p Y+ 


—ieo^ ( or iea^) 
—iea^ ( or iea^) 


W-A^W- 
Y^A^Y- 


iea^ ( or — iecr M ) 
iea^ ( or —iea^) 



Table A. 2 

Z boson vertices 



Vertex 


Factor 


Vertex 


Factor 


Z p v* L v L 

zJ* L i L 


-\ gz{p + p'Y 

^c 2W gz{p + p'Y 


Zfii* R ln 


igzSwip + p'Y 


P P 

Al Xl^/J 


y z {p+p'y 

-h9z{p + p'Y 


P P 

/0* v /0 v 
Al Al 


-haziP + p'Y 
haz{p + p'Y 


P P 

piPtz^ 
pi pi z ix 


\g z o^ (or -\g z o^) 
-\gzC2w^ (or ^gzC2W^) 
yzs 2 w a p (or - \g z s 2 w a^) 


P P ^p, 

~l /— ry 

Pi Pi z^ 

~l /— ry 

P2 P2 Z p 


-{gz^ (or \g Z o^) 
^gzc-iw^ (or -^gzc-iw^) 
-{gzs^o^ (or ^g z s 2 w cr p ) 


w+z p w+ 

Y+Z p Y+ 


-igzcw^ ( or igzcw^) 
-^gzc-iw^ ( or ^gzciw^) 


w-z^w- 

Y^Z^Y- 


igzcy/^ ( or -igzclyd^) 
\gzc-2W^ ( or -^gz^w^) 



A. 4 Mixing in the slepton sector 

As we know, in supersymmetric models, in order to keep the conversation of 
LFV in the lepton sector at tree level, the sources of LFV are assumed to be 
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Table A. 3 



Z 1 boson vertices 



Vertex 


Factor 


Vertex 


Factor 


Z'^l 
Z'Ml 


-| gz>C2w(jP + p'Y 
-%c 2 wgz'(p + p'Y 


Z'pV^R 

z'JVr 


-gz'c z w {p + p'Y 
-igz'Sw{p+p'Y 


P°*P Z', 
Xi x\% 

0* ry) 


-h9z(p + p'Y 

-^9Z'C2W{P+P'Y 

-haz'c^ip+p'Y 


p'°*p'°z', 
x?*x?zi 
x 2 x 2 % 


hgzip+p'Y 
%gz'C2w(p+p'Y 
igz'ctv(p+p'Y 


Wp°z', 
ftptz'v 

pipi 2 '^ 


-\9z>^ (or y Z '^) 

-y Z '^ (or y Z '^) 

y Z '^ (or -\gz^) 


Pi Pi z p. 

~l /— ryl 

P2 P2 n 


{gz*^ (or -y Z '^) 
y z >^ (or -\gz>o^) 
-%gz>C2wv l * (or ^gzC2W^) 


y+z'y+ 

h 1 


y Z 'Ki^ ( or -^gz'Ki^) 


Y^Z'^Y- 


-igz'tit^ ( or ^gz'^l^) 



from the slepton mass terms in the soft-breaking part of the Lagrangian [9,22]. 
For the SUSYE331, there are three mass terms of left-handed slepton, right- 
handed slepton and sneutrios which may independently be sources of LFV. In 
addition, there exists another LFV source original from the Yukawa couplings 
between Higgs and neutrinos. Thus in the SUSYE331 model, there are at least 
four independent sources of LFV and we will parameterize them as follows. 
In each case of supersymmetric particle (sleptons) ip ( ip = II, Ir, ^l,^r), 
we define a corresponding mixing angle 9^ which was defined in [20]. In what 
follows we just remind some general formulas for the review. The mass mixing 
matrices of smuon and stau as well as their sneutrinos can be written in the 
general form of: 



11 2 

' m r m r 



, lib 2 - 



(A.13) 



Mixing angles then can be determined as 



2 

s^ = sin^, c^ = cos^ where s$Cj = — 2 (A.14) 



where = {sl, sr, sp L , sc R } and {rri^ 2 , m^ 3 } are eigenvalues of M 2 ^, accord- 
ing to notations in [20]. In addition, for convenience we denote mj~ instead 
of m^. We always choose m^ 3 < m^ 2 to take the positive values of and 
cr. The mass-eigenstates of sleptons are denoted as {tp 2 , ^3} while the flavor- 
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eigenstates are ip T }. The relation between two bases are: 

= ~ and 4> T = s^ 2 + c^V> 3 - (A. 15) 

B Contribution to r — > fi'j 



Diagrams relating to Cl R are drawn in fig. B.l with no line of Higgs insertion. 

7 (1) 7 (2) % (3) 

w + \w + y + s y + 



\3 \8 



T 



(4) 



7 



7 



Fig. B.l. Dia grams contributing to R . 



Formulas of CI R are: 



Cl 



+ 



(g 2 c L s L ) 

16tt 2 
(g 2 c UL s 



-K^{m\,m 2 



16tt 2 

,2 t a 2 



X 



1 r 

6 



.51 iil x , in? ,m r ,m r ,m? , 

J V A' l L2 i l L2 i l L2 i l L2 / 



nr/ / 2 2 2 2 ~ 2 \ 

-2K 5 (m A , m A , m A , m x ,m VT J 



+ 3m A J 5 (m A , m A , m A , m A , m 

c ur S vr) 



VL2' 



+ 



167T 2 
,2 7 / 2 



6 



-2K 5 (m 2 x , m 2 x , m 2 x , m A , mj^J 



, o 2 t / 2 2 2 2 ~ 2 \ 

+ 3m A J 5 (m A , m x , m x , m x ,m v ) 



16tt 2 162 
(L 2 — > L 3 , R 2 — > R3), 



ov B' i i2 ' l L2 i l L2 i l L2 , 



c 



7 _ g' 2 c R s R 

R 



16tt 2 



/ 2 2 2 2 2 ^ 

°V «' « fl2 ' te' lR2> lR2' 



[R 2 ^ R 3 }- (B.l) 



In the other hand, D 1 gets contributions from diagrams with one line Higgs 
insertion, figs. B.2,B.3 and B.4. There is another class of LFV sources relating 
with neutrino-mediation in which their contributions are very small [21] so we 
will ignore them in our investigation. The D\ R can be separated into three 
parts: 

Dl, R = D^+Dl$ + Df% 
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p'° ^ W p'° >7 ( 2 ) p'° Xb (3) 

\ W+ \ W + \ Y+ < Y+ \ >?a,>?a 



P : (4) f 



?7 

Fig. B.2. Contribution to D^ a) [1-3] and D^ a) [4]. 

where diagrams involving each part are expressed in three figs. B. 2, B.3 and 
B.4. 

For 



D 



7(a) _ 9 Cl-Sl 

16tt 2 



2 

lL2> 



m 



lL2' 



9 c v L s VL 1 

X 



16tt 2 

2 

16tt 2 

,/2, 



2 

1 r 



. 9 c L s L 2 
16vr 2 z " 



2 7 / 2 2 2 2 2 \ 

x t/2 2 2 2 2\ 
4 WriW; ,ffl? ,m r ,f7!r ) 

54 0V ii ' ZlV «L2' «£2' «L2 7 



D 



R 



3 ) 



,/2 



7(a) _ <7 C -R S -R 2 



16tt 2 ^ 2 



t / 2 2 2 2 

« °V ±>> Zoo' «ro' if 



[i? 2 -> i? 3 ] (B.2) 



For D^: 



D 



,7(6) 



16tt 2 



4 t / 2 2 2 2 2\ 

ov A' r^p' f L2 ' y L2' ^£2' 



<? S vr C ur 4 r / 2 2 2 2 2 

-^ H j5(m A , //„, m^, m^, m^ 2/ 



+ 



16tt 2 

2 

5 1 S V L2 C VL2 



x m A /i tan7 [j 5 (m£,m;[,//J,/xJ,m§ M ) 

,t/2 2 2 2 2\,t/2 2 2 2 2 \ 

+ J 5 (m A , m A , m A , // m P ) + J 5 (m A , fi [i p m^) 



+ 



9 S VR2 C VR2 

16vr 2 

,2 2 



x m A /i tan7 J 5 (m A ,m A , //J, /Am? ) 



i T / 2 2 2 2 2\,t/2 2 2 2 2 \ 

+ J 5 (m A , m A , m A , m Pi2 ) + J 5 (m A , \i pl fi p , fi p , m^J 
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p'°i 



'7 



(1) 



VLc 



(4) 



1? jp'° 



VRo 



P 1 



(7) 



'7 



VRo 



(10) 



p'° jP\*%>& 



7j 



(2) 



/ A 1 r c 



p'rpY !if-w+ 



r 



7? tp° 



T 



VRet 



(11) 



Z>0 



Xb 



p'°i 



> 7 



(3) 



P2 \Y + >Y + 



(5) 



^ VRa 

(6) 
1? jp° 

p'i2Pt\w~w + 



/i r c s 



VLc 



(9) 



n'° i 

9 f A B 

p'° \x\,x\ 



7 



(12) 



P°i 



p /0 p° ! X B 



T 



7 



7 

7 (6) 



7 



Fig. B.3. Contribution to L>^° ; [1-10] and D~£ b) [11,12]. 



5- s L c L 2 J 

— mf x - 

16vr 2 3 



7/22 2 2 2 \ 

J 5 {m x ,fi p ,m Lo ,m lL2 ,m lL2 ) 



,2 2 



m A /i p tan 7 J 5 (m A , /i p , , , 
2 



7/2 2 2 
M m BiVpi m 7 



. g' 2 s L c L 2 v 

16tt 2 z " 27 



- m B /i tan7 7 5 (m B , AV m ;~ i2 

- [L 2 ->■ L 3 ], 

2 
9 



2 ^2 _2 2 



lL2 



D ^) = 9^RCR m 2 x _ 

»H2 



7/2 2 2 

Js{m B ,ii ,mr 



2 2 N 

P' "°lR2> m iR2 ,m iR2' 



167T 2 

+ m B /x p tan 7 J 5 (m B , p; 2 ,, m 2 ^ , m} R% , m 2 ^ ) ] 
- [R 2 -> 



(B.3) 



For D^: 



7(c) 



m 



16tt 2 9 



s L c L s R [A T + -fiptmj] + s R c R A„ T + 



„2 „2 



X h 



Bi " L Bi m Bi 



• m L> m U 
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(1) 

Ab Ab 


(2) 


Ab 


(3) 


















(4) 

Ab Ab 


(5) 


Ab 


(6) 




' // 


r y^* 




p° * " ' A \p'° 








(7) 

Ab 


Ab 


(8) 




^ Via ^ ^ ^ 




rI<* /> 




p'°t 7 * 


A 


7^ 




Fig. B.4. Contribution to L»^ (c) 


[1-6] and [7,8]. 





7(c) 



s L c L s R [A r + -/x p tan7] + src^A 



x I 5 (m 2 B ,m 2 B ,m 2 B ,m 



2 m 2 \ 



+ 



slCl c r [A t + -fiptamj] - s R c R A 



x I 5 (m 2 B ,m 2 B ,m 2 B ,m 2 m] ) 



s L c L c R [A r + -/i p tan 7] - s/jCrA 



S L S R A flT 



4- r 2 r 2 4 L 



x I 5 (m B ,m B ,m B ,m~ 



g ^B 

16tt 2 9 



shCb ( s£[A T + ^tan 7 ] + s L c L A^ T ) + c R s 2 L A* T 



x J 5 (m| ,m 2 B ,m 2 B , m}^ , m 2 fl2 } 



s R c R (s 2 L [A T + ^ p tan 7 ] + s L c L A L ^ - s 2 R s 2 L A* T 



x h(m B ,m B ,m B ,m 



lL2 ' m iR3 ' 



+ 



s R c R c L [A T + -/iptan7] - s L c L A 



C R C L^r 



x I 5 (m B ,m B ,m B ,m~ 



m 



lR2 ' 



s R c R (c 2 L [A T + ^/Xptan 7 ] - s L c L A L ^ - s 2 R c 2 L A* T 



x I 5 (m B ,m 2 B ,m 2 B ,m 



lL3 ' ^iRA 1 J ' 



(B.4) 
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C Contributions to Z — » /it 



In this appendix, we draw all the possible diagrams which contribute to the 
effective operator Z — > fir in the limit of assumption given out in [9]. All of 
these diagrams can be applied to the case of Z' boson. 



C.l Contributions to Af R 



Diagrams contributing to A Z L R are shown in fig. C.l. The formulas are: 

,0 /o /0 /0 

' ' 1 ' (2) P : P : (3) ' : (4) 
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Fig. C.l. Diagrams contributing to ( or ^) (first, second and third rows) 



i Z R {a) ( or A R Z ' iah 
implies i,j = {B, 3, 8} and i / j. 



and ( or ^ w ) (fourth row). Here we denote JETjg G while A 



,Z(a) 

1 L ~~ 



n 2 r 2 1 
x^x (1 +C27 ) 



167T 2 4 

2 t / 2 2 2 2 2 \ ot/2 2 2 2^ 



+ \ S V R C V R , 



n 2 r 2 1 



X 



2t/2 2 2 2 2 \ ot/2 2 2 2\ 

-V P M m \, ™ x , ti p , ti p , m Pm ) - 2J 4 (m x , m x , \i p , m ) 
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, . g 2 cL 11 



2 7 / 2 2 2 2 2 \ ot/2 2 2 2 \ 

x -Li p J 5 {m x , m x , Li p , Li p , m~ lL2 ) - 2 J 4 (m A , m x , fi p , m^J 

, 2 ( j I 2 2 2 2 \ 2r / 2 2 2 2 2 \\ 
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X 16?r 2 X 4 i i 



C2o 



X 



2r / 2 2 22 2 \ 7/2 2 2 2 \ 
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/ 2 t / 2 2 2 2 2 1 / 2 ' 2 

[Li p h[m B , m B , ii p , Li p , mj L2 ) - h[m B , m B , Li p , 

(1/2 — > -^3) -^2 - > -^3)- 



(C.l) 



Ft 



, , 9' 2 t 2 c 2 w 2 

X 



x {vpMmli m Bi tfpi m l R2 ) + 2M m Bi m Bi V% rnf R2 ) 

2 [ 2i / 2 2 2 2 2 \ 

- ™b [Vph{m B , m B , Lip, Lip, m- lR2 ) 

- h{m 2 B ,m 2 B , Li 2 p ,m 2 R2 )\] - [R 2 ->■ i? 3 ] 



(C.2) 



For A 



Z(b,c) 



Formulas for 



A Z L ib) = (s L c L ) 



m r C W h f 2 2N 



2|T/2 2 2 2 2 \ 

S R ^5("l A > ^> m I L2 1 m i R2 ) 
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(!) 

P KA 3 ,A 8 ) 



(2) 



(Ab, A 8 ) 
A 3 



(3) 



o A 



,o r 



Fig. C.2. Diagrams contributing to ^ (left side) and -A^jj (right side) in 
SUSYE331. 
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16tt 2 V 2 27 v 7 ^ p 

7/222 2 2 \\ 

+ J 5 (m B ,^m fM ,m fiB ,m r JJ 
+ J 5 (m 2 B ,fi 2 p ,m 2 L2 ,m 2 m ,m 2 lR3 ))] - (L 2 ^ L 



( T I 2 22 2 2 

{J 5 {m B ,fi pJ n p ,m !L2 ,m lR 



( t I 2 2 2 2 2\ 

[J 5 {m B ,^ p ^ p ,m iL2} m !mJ 



(C.3) 



where m r = F T x t>'/ a/2 is mass of the tau and V = w wea k = Vv 2 + v' 2 in the 
SUSYE331 model. We also have formula of A B ^: 



A z(b) , s mlt 2 ^ 1 2 2 



16tt 2 V 2 9 



X 



-s £ ( J 5 {m B ,n p ,fi p ,m iL ^m iR2 



i 7 / 2 2 2 

2/7/2 2 
c L \J 5 [,m B ,n p ,fi p ,m 

+ J 5 (m 2 B ,fi 2 ,m 



«2' 



2 2 
Z~L, ' '"fin 



2 2 



(i?2 -> #3) 



(C.4) 



Formulas for A z ^: 
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Z(c) 
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X ^ s 



27/2 2 2 2 2" 
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Il 3 ' II 



,27/2 2 2 2 2 \ 

^ i_R 3 ' 'h 3 ' Zi 3 ' «Z, 3 7 
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- 1 4 - 4 



2 t / 2 2 2 2 2\ 

s R J 5 (m A , m f , m- , m f , m f ) 



2 t / 2 2 2 2 2 \ 1 



16tt 2 I/ 2 108 



< 4 



2 t / 2 2 2 2 2 \ 



«2' 



, 2 t / 2 2 2 2 2 \ 

+ c fl J 5 (m B ,m ffls ,m fiiQ ,m rio ,m fio ) 



- c r 



2 t I 2 2 2 2 2\ 

/ 2 2 \ 2 t / 2 2 2 2 2 \ 

" [H ~ c L J [s R J 5 (m B , m lR2 , , , ) 

2t/2 2 2 2 2\i 

+ c R J 5 (m B , mj R3 , , , ) j j 



(C.5) 



A Z R ic) = (s R c R ) 



x 



^t 2 ^ J_ 2 2 

16tT 2 \/ 2 12 W 
2t/2 2 2 2 2 \ 

s T J 5 (m R ,mr ,m? ,m r ,m r 



+ c 2 L J 5 (m B ,m 



i i i \ 
, ,m? ,m r , m,~ ) 

l R 2 l R 2 l L 3 1 <L 3 ' 



+ 4 



2t/2 2 2 2 2\ 

s T J 5 (m R ,m r ,m T ,m r ,m r ) 



, 2 t / 2 2 2 2 2 \ 

+ c L J 5 (m B , , m- lR3 , m fj . s , ) 



2t/2 2 2 2 2 \ 

s L J 5 {m B , m ; - K _ , m L , m f , m f ) 



1 Ir 2 



2 



+ (4 - 4) 

2t/2 2 2 2 2 \ I 



(C.6) 



C.2 Contributions to C^ R 



For Cf in fig.C.3. The formulas for these two quantities are written as the 
below: 



g 2 1 

CE = (c UL s UL ) x — x — 
.2 x 



c / 2 2 2 2 2 



+ ( ° lSl) X 16V 2 X I8° 2W ~^ 5 ( m '' m L' m L' m L' m D 
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w+ s w + 



(1) 



(2) 



r 



(A3, As) ( 3 ) 



(z') 

V (4) 



>Z{Z>) 
(5) 



(6) 



Y+ 



T 



VRa 



VRa „>' 



Fig. C.3. Diagrams contributing to Cf R {Cf R ). Only the last gives contribution to 
C R (Cf)- The first dia 

gram only contributes to while the fifth only contributes 

<z' 

L 



to C?'\ 



( \ g2 1 2 

+ \Cu L 8u L ) X X Y2° W 

, o 2 7 / 2 2 2 2-2 

+ 3m x J 5 {m x ,m x ,m x ,m x ,m 



-2K 5 (m 2 x , m 2 x , m 2 x , m 2 x , m 2 ) 



VL2' 



T ( Cj/p St, 



C2W 



' VR VR > 16tt2 12 

io2t/2 2 2 2 - 

+ 3m x J 5 {m x , m x , m x , m x , m 



2K 5 (m 2 x ,m 2 x ,m 2 x ,m 2 x ,m 



VR2< 



g' 2 1 
16tt 2 324 
[L 2 — >■ £3, -R2 — >■ -R3] , 



+ (c L s L ) X 



5 w0 



7>- / 2 2 2 2 2 \ 

A5 m R ,mr ,m ; ,m? , mr 

°v HI l L2 i l L2 i i L2 > l L2 i 



(C.7) 



c r = (crSr) x 

- [R 2 ->• #3 



5' 2 1 2 

16^ X 36^ 



K^(m 2 n ,m 2 , mf ,mf , m 



2 * 



(C. 



We note that because Z boson couples much weakly to right-handed neutrinos 
so the diagram 5 in the fig.C.3 give suppressed contribution to . In contrast, 
the case of Z' boson is different, it weakly couples with W ± but non-negligible 
to right-handed neutrinos. So for the C z , we neglect the first diagram and 
keep the fifth. This conclusion is held in the case of D z and D z . 



C.3 Contributions to Df R 



For Df R , we have Df R = D^ R + D^} R . They are presented by diagrams in 
figs. C.4 and C.5. 
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(1) 
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Z?p° 

7$-\ 7& 



(2) 



P°i 



(3) 



Pi Pi W~i w+ pT>Pi\W-W + Pi pi \W~W + 



fZ')(4) (Z) ( 5 ) 
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( Z ') (8 
^P°f(A B ,A 8 ) 



Z' 



P P 



p'°2 p°\ 



p'° S p° 



(10) 



T C ^ It '' V 



T 



'Z(Z') 



>Z(Z>) 



Fig. C.4. Diagrams contributing to D^ b ^ (D^ ^) (two first lines) and D^ b ^ {D^ ^) 
(the last line). Noting that the first diagram only contributes to while the 

sixth only contributes to ^ . 



Formulas for D z ^: 
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, 9 2 1 
(slCl)y^ x g x /i p m A tari7 

ot/2 2 2 2 2\,t/2 2 2 2 2 \ 

2 J 5 (m A , // // // mr J + J 5 (m A , m A , // n m,~ ) 



+ (slCl) x -yU P ^Atan7 



x c 2W 



2 r / 2 2 2 2 2 \ 
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9' 



,,- !7j— 7 x ^ (/i p m A tan7) x 



x c 2 vf 

+ ( S U L C U L 



o 7/2 2 2 2 2X17/2 2 2 2 2- 

2J 5 (m x ,m x ,m x , fi p ,m 0ii3 ) + J 5 {m x , m x , [i p , fi p , m 5fi2 , 



1— x (/i p m A tan 7) x - 
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X 



(-1 + 2s 2 w ) [2 J 5 {m 2 x , /i 2 p , fi 2 p , fi 2 p , m 2 L2 ) + J 5 {m 2 x , m\, /i 2 p , fi 2 p , m 2 



2 - 



. . 2 To T ( 2 2 2 2 2 \ , t / 2 2 2 2 2\ 

x s w [2J 5 {m x , n p , n p , n p , m Vm ) + J 5 {m x , m x , fi p , fi p , m Pja ) 

g' 2 1 
+ (s L c L )^-j x 27 x /V 7l £ tan 7 

2 r / 2 2 2 2 2 \ 

x c 2 wm lL2 h{m B , n p , mj M , m^m-^ 
+ 16t^ X 54 X /V 7l - Btan T 

o 7/2 2 2 2 2 \ 1 T / 2 2 2 2 2- 

2J 5 (m B , fi p , n p , fi p , m~ lL2 ) + J 5 (m B , m B , /x p , /i p , m^, 
— (L 2 — >■ L 3 , i? 2 -» i? 3 ), 
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Z(6) 



5 ' 2 1 

-(sflCiijjg-j x —m B fi p taii7 



'-4s 



2 \ 2 T I 2 

w) m i R2 h{m B , fi 



2 2 



+ 2J 5 (m 2 B , fi 2 p , /I 2 , /i 
-(R 2 ^R 3 ). 



2 2 N 



t 1 2 2 2 



"pi m i 



(CIO) 



Ab Aff 

^z i P ° //'! z 

(z') ' (z') 

Ab A_b 



^z V p or 2 

(z') ' (z') 

Fig. C.5. Dia grams contributing to iJ-^Lli^ 
Formulas for D^ R : 
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+ 4 J 5K, m B> m L' m L- m L) )1 - ( L 2 -> ^ 3 ), (en) 



Dr = -{srCr) x -^m B p p tan-f x — 

X 



(1 - 24,) (4 J 5 (m|, m\,m} L2 ,m} L2 ,m} R2 ) 

,2t/2 2 2 2 2\\ 

+ c L Mm B ,m B ,m iL3 ,m lL3 ,m !R2 ) ) 

, o 2 / 2 t / 2 2 2 2 2 \ 

+ 2s w (s L J 5 (™b, m B , m lL2 ,m !R2 ,m iR2 ) 

+ clJ 5 (ml,ml,ml 3 ,m] R2 ,m 2 !R2 ) )] - (i? 2 -)■ i2 3 ). (C.12) 



D Contributions to Z 7 — > pr 



£>. 1 Contributions to A 1Z ' R 



To determine the values of R , A XZ R and Af^, we use techniques mentioned 
in [9] . From formulas of covariant derivatives of neutral Higgs in the appendix 
A. 2, it is easy to see that two terms relating with Z and Z' bosons appearing 
in these covariant derivatives are different from each others one factor (— 1). 
For A\ Z R which relates with p° and p'° we have A IZ R = A^' R {a) + A Z L J ] + A Z L '§ 
and ^4f( fl ) = ( m z/ m z')A 1 I ^ R -j + A^' R y This leads to the results: 

A lZ'{a) _ A Z(a) 
^L,R — L,R 1 

A Z'(b) _ _ A Z{b) 
L,R — ^±L,R ) 

<J? = -A Z L % (D.l) 



D.2 Contributions to A 2Z R 



Diagrams contributing to A 2Z ' R are quite similar to those shown in flg.C.l. 
There is a interesting point in the SUSYE331 model that both Higgs x an d 
x' do not couple to leptons and sleptons. As a consequence, A L R give contri- 
butions from only the class of diagrams depicted in fig.D.l where p°, p'° and 
Z boson in flg.C.l are correspondingly replaced with x%i X2 an d boson. We 
use the equivalence role between p° <r-> x%. an d p'° -H- xf f° r the calculation ( 
see App.A.2). The results are: 
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• is » ■ « ■ ► , « 

T V ^ X /* 
- > — ' 




' c 







(2Z f ) (2Z') 

Fig. D.l. Diagrams contributing to ^4^ (first, and second rows ) and AY ' (third 
row). Here we denote H® G {xlL X?} while Ay implies i,j = {B, 8} and i ^ j. 
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(D.3) 



Contributions to C^' R 



Diagrams contributing to C^' R are those from 2-6 in the fig. C.3. Comparing 
with the case of the Z boson we easily deduce the formulas as: 



Cf — — {c Vl s Ul ) x 



9' 



1 

X —C 2 W 
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(D.4) 
(D.5) 



D.^ Contributions to Df' R 



Contribution to Df R can be deduced from diagrams shown for R in figs. 
C.4 and C.5. We also write D? R = £>f { R } + £>f { R c) . From fig. C.4 we can deduce 
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formulas to determine D z '^: 
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E Contribution from B^ L R R to r — > 3fj, 



Contributions to B^ L R R arise from the diagrams in fig. E.l. The formulas are: 
r t>L a fi t VR a V> t lh a [It h a 
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>W+ W + 
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Xb Xb 



Xb Xb 
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Xi Xi 

— — 

Xb Xb 



Xb Xb 



Fig. E.l. Diagrams contributing to B^ L ' R (first and second rows) and B^ R ' R (third 
row). Xi and Xj (i 7^ j in each above diagram) are gauginos in which Aj and 
Xj G {Xb, A3, As}. 
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